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Abstract. DFT calculations have been used to provide insights into the origin 
of the colossal positive and negative thermal expansion in Ag3[Co(CN)e]. The 
results confirm that the positive expansion within the trigonal basal plane and the 
negative expansion in the orthogonal direction are coupled due to the existence 
of a network defined by nearly-rigid bonds within the chains of Co— C— N— Ag-N- 
C-Co linkages. The origin of the colossal values of the coefficients of thermal 
expansion arise from an extremely shallow energy surface that allows a flexing 
of the structure with small energy cost. The thermal expansion can be achieved 
with a modest value of the overall Gruneisen parameter. The energy surface is 
so shallow that we need to incorporate a small empirical dispersive interaction 
to give ground-state lattice parameters that match experimental values at low 
temperature. We compare the results with DFT calculations on two isostructural 
systems: H3[Co(CN)g], which is known to have much smaller values of the 
coefficients of thermal expansion, and Au3[Co(CN)g], which has not yet been 
synthesised but which is predicted by our calculations to be another candidate 
material for showing colossal positive and negative thermal expansion. 



PACS numbers: 61.43.Fs, 62.20.Dc, 63.50.+X 

1. Introduction 

Recently we reported diffraction measurements showing colossal thermal expansion 
in the network structure of Ag 3 [Co(CN) 6 ] [JJ. Within that paper we presented data 
showing a positive linear coefficient of expansion within the basal plane of the trigonal 
unit cell of a.\ — +132 MK _1 , and a negative coefficient for the orthogonal c axis of 
a 3 = -130 MK" 1 (1 MK _1 is equivalent to 1 x 1CT 6 K" 1 ). To highlight just how large 
these values are, we compare values for a of number common materials in table [T] In 
the recent paper pQ we briefly summarised some of results from an ab initio Density 
Functional Theory (DFT) study; in the present paper we report the results from this 
study in more detail, together with an extension to the related materials H3[Co(CN)e] 
and Au 3 [Co(CN) 6 ]. 

There are two main issues involved in understanding the thermal expansion of 
Ag3[Co(CN)e]. First is to understand whether the positive expansion in the trigonal 
basal plane and the negative expansion in the orthogonal direction have a causal 
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Table 1. Coefficients of the linear thermal expansion tensor, a n , for a number 
of common materials. In none of these examples is the thermal expansion caused 
in any part by the existence of a phase transition. Units of a are MK . 



Material 


Temperature (K) 


ai «2 


013 


references 


Ag 3 Co(CN) 6 


10-500 


132 


-130 


m 


d-H 3 Co(CN) 6 


4-300 


14.8 


-2.4 


m 


ZrW 2 O g 


0.4-430 


-9.1 




mm 




430-950 


-4.9 




mm 


low cordicritc 


600-1050 


6 5 


-0.6 


m 


Cd(CN) 2 


150-375 


-20.4 




6 


/3-quartz 


900 


1.9 


-1.1 


.7 


NaCl 


293 


39.6 




A 




Figure 1. The crystal structure of Ag3[Co(CN)s], space group P31m, showing 
one unit cell viewed in perspective (left), and as viewed down the crystallographic 
[001] axis (right). The atoms are labelled in the diagram. We note that the Ag 
atoms are the large (red) spheres at z = 1/2, the Co atoms are the small spheres 
at the origin of the unit cell (dark blue), the N atoms are the spheres closest to 
the Ag atoms (light blue), and the C atoms are the spheres closest to the Co 
atoms (black). Points to note are the existence of CoCg octahedra located on the 
corners of the unit cell, and the Co— C— N— Ag-N— C— Co chains lying across the set 
of {100} planes. Colour is only available in the on-line version. 

relationship. An inspection of the crystal structure, figure [l] suggests that the origin 
of the coupling between the positive and negative expansions lies in the way that 
the structure forms networks consisting of chains of Co-C-N-Ag-N-C-Co linkages, 
connecting CoCg octahedra via N Ag N linear groups. If the Co... Co chains remain 
nearly-rigid and nearly-linear, it is an inevitable consequence of the existence of this 
network of chains that any expansion in the basal plane has to drive a contraction 
of nearly equal magnitude in the orthogonal direction. However, this fact needs to 
be quantified. The second issue concerns what features of the material give rise to 
colossal values of the coefficients of thermal expansion. 

Issues such as these are within the domain of atomistic simulation. In the first 
instance, we judged that empirical models for Ag 3 [Co(CN) 6 ] would probably be hard 
to develop (for reasons that include the paucity of appropriate experimental data for 
tuning parameters, and the fact that the bonding is not likely to be described by 
the simple standard functions that would provide the necessary starting point), and 
thus we have started directly with DFT methods. These are described in the next 
section, and the main results are described in the subsequent section. We also report 
a set of calculations for the isostructural material H3[Co(CN)e], which has the same 
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coupled positive and negative expansion but with coefficients that are closer to most 
other materials (table[T]), and for Au3[Co(CN) 6 ], concerning which there are currently 
no experimental data but which is predicted from our calculations to have similar 
behaviour to Ag3[Co(CN)g]. The final section reports a simple analysis based on the 
calculations, showing that the colossal thermal expansion in Ag3[Co(CN)g] can be 
accounted for with a modest value of the overall Griineisen parameter. 

We would like to point the reader to a parallel paper that reports the results of a 
neutron total scattering study of Ag 3 [Co(CN) 6 ] [5]. It contains additional discussion 
of the crystal structure, and presents a number of results that support the work of the 
present paper. 

2. Methods 

2.1. DFT calculations 

The DFT calculations were performed using the academic version of the CASTEP 
code, v4.2 |10j operating in energy minimisation mode. This approach uses a plane- 
wave representation of the electronic basis states, and represents the inner atomic 
electronic orbitals through the pseudopotential approximation. We used both norm- 
conserving (NCP) and ultrasoft (USP) pseuopotentials in various calculations, always 
using the generalised gradient approximation (GGA) to the exchange-correlation 
energy. The NCP pseudopotentials were generated using the OPIUM code [TTj.and 
the USP pseudopotentials were used as supplied with CASTEP with the exception 
of Co, which was kindly generated for us by Dr Victor Milman (Accelrys). In all 
our calculations we relaxed the crystal structure to give the lowest energy; this 
always involved relaxing the atomic coordinates, and in various calculations we cither 
performed a full relaxation of the lattice parameters, or a partial relaxation, or we 
held the lattice parameters at fixed values. 

Anticipating the results we present below, we found that DFT is unable to 
reproduce the correct ground-state structure of Ag3[Co(CN)g] for reasons that we 
subsequently realised point to some scientific insights. Thus we tried a number of 
approaches, and in this paper we report the calculations on Ag3[Co(CN)g] using both 
the NCP and USP models; the calculations on H 3 [Co(CN) 6 ] and Au 3 [Co(CN) 6 ] were 
performed using the NCP model only. The plane- wave energy cutoffs used in the NCP 
and USP calculations were 680 and 470 eV respectively, obtained as good values from 
convergence tests. As is usual for this type of calculation, the Bloch functions were 
sampled over a grid of wave vectors using the standard Monkhorst-Pack grid [T2] , in 
both cases finding that converged results were obtained using a 3 x 3 x 3 grid. 

The trial crystal structure used in each simulation was that determined by Pauling 
and Pauling (T3] for Ag 3 [Co(CN) 6 ], as confirmed by our recent neutron diffraction 
measurements [1 . Where we adjusted lattice parameters (e.g. for running calculations 
across a grid of pre-set lattice parameters) we used the same set of experimental 
fractional atomic coordinates for the starting structure. The trial crystal structure 
had space group P31m (see table [2]). There is no automatic method for constraining 
the crystal symmetry during the energy relaxation process within CASTEP, but the 
output results showed that the symmetry did not change through any of our relaxation 
calculations. The simulations on H3[Co(CN)6] were started with the experimental 
lattice parameters and fractional atomic coordinates. In the absence of experimental 
data, the simulations on Au3[Co(CN)g] used the trial structure of Ag3[Co(CN)e] as 
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Table 2. Results for the calculated ground state of Ag3[Co(CN)e] from the 
different models using the CASTEP calculations: norm-conserving pseupotentials 
(NCP), ultrasoft pscudopotcntials (USP), and compared with experimental data 
fo temperatures of 10 K [9] and 300 K jT|. Note that the asymmetric atomic basis 
is defined with coordinates Co (0,0,0), Ag (1/2,0,1/2), C (x,0, z), N (x,o, z), 
with values for the atomic fractional coordinates x and z given here. 



Experiment (10 K) Experiment (300 K) NCP USP 



a (A) 


6.754 


7.031 


7.764 


7.697 


c(A) 


7.381 


7.117 


6.485 


6.681 


C x 


0.222 


0.210 


0.199 


0.198 


C z 


0.156 


0.158 


0.174 


0.167 


N x 


0.339 


0.335 


0.317 


0.317 


N z 


0.264 


0.269 


0.286 


0.278 


C-N (A) 


1.122 


1.182 


1.205 


1.179 


Ag-N (A) 


2.053 


2.012 


2.062 


2.045 


Co-c (A) 


1.890 


1.856 


1.916 


1.889 


Ag-Ag (A) 


3.377 


3.516 


3.882 


3.849 



the starting point. 

2.2. Computational details 

The CASTEP jobs were all run on the high-throughput computing grid of the 
University of Cambridge, CamGrid [TU [TS]. This consists of a large number of 4- 
core (dual processors, with each processor having dual cores, and with 2 GB memory 
per core) shared-memory machines. Jobs were scheduled using the Condor middleware 
[T6] . Each job was restricted to a single node (one node being a 4-core machine), and 
run as a 4-processor parallel task. CASTEP was compiled with an Intel F95 compiler 
(v9.1), and built for the LAM flavour of MPI [17 . By selecting the coll.shmem 
SSI module for LAM, all inter-process communications used shared memory. This 
avoided the computation making unnecessary use of the network stack, which would 
have degraded communication speed. 

3. Results 

3.1. Ground state calculations on Ag^[Co(CN)(,] 

The results of complete structure-relaxation calculations for both USP and NCP 
models are presented in table [2j The striking point is that there is a substantial error 
in the values of the a and c lattice parameters - up to 15% - although the Co-N, Ag-N, 
and C-N bond lengths are in good agreement with experiment, with errors of around 
2% (maximum error is for the Co-C distance in the NCP calculation and is 3.6%). 
The discrepancies in individual bond lengths are typical for this type of calculation 
(the community experience is that GGA usually leads to a small over-estimate of 
bond lengths), but the large errors in calculated lattice parameters are not expected 
for the methods we used. We believe that the large discrepancy between experiment 
and calculations is actually providing some critical information about the physics and 
chemistry of Ag 3 [Co(CN) 6 ], suggesting that an important feature is missing from the 
DFT calculations. Anticipating results described later in this paper, we note that our 
calculations on H 3 [Co(CN) 6 ] gave results that are much closer to experiments (within 
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Table 3. Mulliken atomic charges (Q) and Mulliken bond orders (BO) 
calculated for Ag3[Co(CN)g] (first four columns of data) and H3[Co(CN)6] and 
Au3[Co(CN)6] (last two columns). NCP and USP correspond to the use of norm- 
conserving pseupotentials and ultrasoft pseudopotentials respectively, and MIN 
and EXP refer to the lattice parameters of Ag3[Co(CN)g] after minimisation and 
with the experimental values respectively. 





NCP_MIN 


USP_MIN 


NCP_EXP 


USPJEXP 


H 3 [Co(CN) 6 ] 


Au 3 [Co(CN) 6 ] 


Qc„ 


1.13 


1.15 


1.16 


1.19 


1.32 


1.18 


QAg/H/Au 


0.68 


0.61 


0.65 


0.58 


0.36 


0.54 


Qc 


-0.01 


-0.02 


-0.01 


-0.02 


0.05 


0.00 


Qn 


-0.52 


-0.48 


-0.51 


-0.47 


-0.45 


-0.47 


BO Co-C 


0.25 


0.20 


0.25 


0.19 


0.22 


0.25 


BO C-N 


1.80 


1.84 


1.80 


1.84 


1.77 


1.81 


BO Ag/H/Au-N 


0.33 


0.35 


0.34 


0.37 


0.43 


0.41 


BO Ag-Ag/H-H/Au-Au 


-0.01 


-0.03 


0.00 


-0.02 


0.00 


0.00 



the usual GGA level of agreement), suggesting that the discrepancies are due to some 
specific feature of the Ag cations. 

Given that the DFT calculations are giving reasonable results for the individual 
bond lengths, we repeated the calculations in which the lattice parameters are held 
at their experimental values and atomic coordinates were allowed to relax. For the 
NCP and USP models, the energy differences are 0.111 and 0.169 eV per formula unit 
respectively (0.037 and 0.056 eV per Ag atom respectively). This is actually quite 
a small difference, and suggests that whatever is missing from the DFT calculations 
is not large, but nevertheless appears to be significant in the present case because 
overall energy differences are small. It is also interesting to note that the DFT 
results do not show any appreciable rearrangement of the atomic charges or change 
in the bonding characteristics when changing the lattice parameters — the calculated 
Mulliken charge distribution and Mulliken bond orders obtained directly from the 
CASTEP calculations [HJ [19] are given in table [3] and show no significant changes 
with large changes in lattice parameters. 

Our hypothesis is that the large discrepancies between the experiment and 
calculated ground state structures can be traced to the Ag...Ag nearest-neighbour 
(non-bonded) distance, which is half the value of the a lattice parameter. This is 
longer in the calculation than experiment by 0.5 A. The Ag...Ag distance at low 
temperature - 3.38 A - is of the same size as twice the Van der Waal's radius (1.72 
A), and thus might be considered to be rather short in view of the fact that this 
is a charged cation (DFT charge value is of order 0.6 positive electron units). The 
electron density obtained from the NCP calculations is presented in figure [2] This 
shows that there is no covalent bonding between the neighbouring Ag cations, although 
covalent bonding is seen along the Co-C-N-Ag chains. Thus we believe that the the 
expansion in the basal plane in the calculation relative to the experimental structure 
can be understood as arising from a repulsion between the neighbouring Ag atoms, or 
cquivalently, due to the lack of an additional attractive interaction in the model. In 
short, there is some aspect of the Ag...Ag interactions that is not properly reflected 
in the DFT calculations, albeit one that will not give large changes in energy (recall 
the comparison between the calculations with complete structure relaxation and those 
held at the experimental lattice parameters above). 

At this point we note the work of Pyykko and co-workers on metallophilic 
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Figure 2. Electron density contour plots for two planes in the crystal structure 
of Ag3[Co(CN)6]. The vertical plane (llO) contains the Co-C-N-Ag linkages 
and shows the strong covalent bonding within the C— N group and additional 
covalent Co— C and Ag-N bonds as electron density between the bonded atoms. 
The horizontal plane ((001) plane at z = 1/2) contains the nearest-neighbour Ag 
cations, and shows no electron density between these atoms. The labelling of the 
atoms is as given in figure [l] Colour is only available in the on-line version. 

interactions involving group 11 elements [20j EH E2]- These appear to arise due 
to electron correlation effects that are not included within the DFT approach [23 , 
and can be represented as an enhanced dispersive interaction. Based on this, we 
have attempted to correct for the discrepancies between calculation and experiment 
in Ag3[Co(CN)g] by adding a post hoc dispersive interaction between the Ag atoms 
of the form — CY~ 6 . We find that the minimum energy state with lattice parameters 
closest to the experimental low-temperature values can be achieved using a value of 
C of 80 eV A 6 per atom pair. We should remark that in general post hoc corrections 
would not be expected to give consistent results, because we minimise the atomic 
coordinates without the correction. However, in the present case the calculations 
show that the bond lengths remain more-or-less constant independent of the values 
of the a and c lattice parameters, and thus we do not expect the inclusion of the 
post hoc correction to significantly affect bond lengths. Of course, we could have 
also included dispersive interactions involving the C and N atoms, but based on our 
calculations of H3[Co(CN)6] presented below we have compelling evidence that the 
Ag...Ag interactions are the critical ones here. In any case we anticipate that any 
effects of the C and N atoms will be incorporated into the Ag...Ag interaction we 
included through the value of C we used. Thus we conclude that Ag...Ag interactions 
that are not captured within the DFT method are required to obtain the correct 
ground state of Ag3[Co(CN)g]. However, these are not large and are only noticeable 
because of the high degree of flexibility of the network that defines the crystal 
structure. We show below that this additional interaction actually softens the energy 



DFT study of colossal thermal expansion in Ag 3 [Co(CN) 6 ] 7 




6.0 6.5 7.0 7.5 8.0 8.5 6.0 6.5 7.0 7.5 8.0 8.5 

a lattice parameter (A) 



Figure 3. Contour plots for the grid of calculations for Ag3[Co(CN)5] holding 
a and c at fixed values, obtained using with norm-conserving pscudopotentials 
without (top left) and with (top right) the post hoc dispersive interaction. The 
line of points (red in on-line version) show the locus of experimental data. The 
lower curves trace the energy of the locus of points along the bottom of the valley 
in the contour plots. 

surface, and is implicated in the ability of this material to show colossal positive and 
negative thermal expansion. 

3.2. Energy surface for Ag 3 fCo(CN) 6 ] 

In figure [3] we plot the energy surface computed for a grid of fixed values of a and c for 
the model with norm-conserving pseudopotentials, allowing relaxation of the atomic 
coordinates for each calculation. The results are shown with and without the post hoc 
dispersive interaction. The key result is that the energy surface consists of a very steep 
valley with an extremely shallow 'floor'. The post hoc additional dispersive interaction 
has little effect on the overall shape of the valley, but shifts the minimum point along 
the floor. Figure [3] also contains the locus of experimental data for a and c for various 
temperatures [T]; the position of this locus is very close to the calculated valley floor, 
confirming that the DFT calculations are capturing the correlation between the values 
of a and c seen in experiment. In figure [3] we also plot the energy from from the locus 
of minima points obtained by holding fixed values of a and b and allowing c to relax. 
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Figure 4. Plots of the angles (top) and orientations with respect to the (001) 
plane (bottom) of the nearest-neighbour bonds in Ag3[Co(CN)s], calculated using 
fixed values of a and allowing c and the atomic coordinates to relax. The angles 
are defined such that the angles in the cartoon, which corresponds to the structure 
for a = 6.0 A, are both acute, and the orientations are given with respect to the 
horizontal line (see cartoon). 



3.3. Structure analysis 

We note here that the Co-N, C-N and N-Ag bond lengths do not change significantly 
across the range of the calculations. Maximum changes in bond lengths for the Co-N 
and N-Ag bonds are 1% and 0.8% respectively for the range of calculations in which 
we held the a lattice parameter fixed at values between 6 and 8.5 A, and allowed 
relaxation of the c lattice parameter and of the atomic coordinates. The change in 
the C-N bond length across the range of calculations is even smaller, 0.1%. 

Thus the largest changes in the structure for large changes in the lattice 
parameters are in the orientations of the bonds, or equivalently in the bond angles. 
In figure [4] we plot the calculated Co-C-N and C-N-Ag bond angles (top), and the 
calculated orientations of the Co-C, C-N and N-Ag bonds with respect to the (001) 
plane ( bottom) , all as a function of the lattice parameter a obtained with relaxation of 
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Figure 5. Plots of the DFT band structure of Ag3[Co(CN)s] computed at the 
values of lattice parameter a = 7.7652 A of the minimum DFT energy structure 
(left) and for a = 6.0 A (right; lattice parameter c and atomic coodinates were at 
relaxed values). 

c and the atomic coordinates. Note that a linear Co-C-N-Ag chain will correspond 
to the three orientations being equal. Thus we note that for large values of the a 
lattice parameter we have a greater tendency towards forming linear chains, but on 
decreasing the value of a we see significant crumpling of these chains. The orientation 
of the Co-C bond is actually determined by the shape of the CoCg octahedron, and 
the fact that this changes the least - by no more than 5° - is consistent with the 
fact that the CoCg octahedra retain their basic shape across the range of calculations. 
The C-Co-C angle (not shown in figure [4| varies in a more-or-less linear manner from 
84° at a = 6 A to 92° at a = 8.5 A. In order to accommodate the large decrease in 
the value of a, and the concomitant increase in the value of c, it is essential that the 
chain be aligned closer towards the vertical direction. Thus there is a clear need for 
the chain to buckle, as seen in figure [4] It is interesting to note that the two chain 
bond angles change by similar amounts (figure [i] to p), so that the N-Ag bond changes 
its orientation by the greatest amount (figure |4j bottom), rather than the buckling 
being accommodated more by bending of one bond angle than the other (note that 
the N-Ag-N linkage is constrained to be linear). 
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3.4- DFT Band structure 

The DFT band structure of Ag3[Co(CN)g] is showin figure [fj] for two different values 
of lattice parameter a, namely 6.0 A and 7.7652 A (the lattice parameter of the DFT 
minimum energy structure), computed on structures for which the lattice parameter 
c and the atomic coordinates had been relaxed. Two points should be noted from 
the data. First, the DFT band gap - the existence of which is consistent with the 
experimental observation that Ag 3 [Co(CN) 6 ] is an insulator - increased from 2.84 eV 
to 3.37 eV on increasing a from 6.0 A to the value of the DFT energy minimum 
structure. The band gap increases further to 3.70 eV on increasing a to 8.4 A. 
Although DFT is known to significantly underestimate the true band gap, trends in the 
behaviour of the band gap may be considered to be qualitatively reliable. Thus, whilst 
the changes in the DFT band gap are not large, there is nevertheless an appreciable 
lowering of the band gap as the Ag atoms come closer together, and as the C...C 
chains become more vertically inclined. The second point to note is that the splitting 
of the top two valence bands increases significantly on decreasing a, with one band 
developing a greater variation with wave vector. These two bands are degenerate at 
the r and A (k = 0,0, 1/2) points (and along the T-A line), and there is maximum 
splitting along the L-M and H-K lines (L, M and H are k = 1/2,0,0, 1/2,0,0 and 
1/3, 1/3, 1/2 respectively [M]). 

3.5. Aside: orientation of the CN molecular ion 

From a purely electrostatic viewpoint, the Co cations might be expected to have 
N rather than C as its neighbouring atom, particularly as the N component of the 
CN~ molecular ion contains most of the charge (as seen in table pjj. Indeed, all 
known structural analogues, such as La[Ag(CN)2]3 and KCo[Ag(CN)2]3, have the CN 
molecular ion rotated by 180° relative to Ag3[Co(CN)e]. This raises the question of 
whether the orientation of the CN molecular ion in Ag3[Co(CN)g] is in a genuine 
cqulibrium arrangement, or whether this structure is the metastable result of the 
growth process. We therefore calculated the energy for the case where we swapped 
the positions of the C and N atoms, giving Co-N-C-Ag-C-N-Co chains instead of the 
Co-C-N-Ag-N-C-Co chains found in the experimental structure. Our calculations 
showed that the the experimental structure has the lower energy difference, with an 
energy difference of 1.03 eV per unit cell from the DFT calculations with the norm- 
conserving pseudopotentials. 

3.6. Calculations on H 3 [Co(CN) 6 ] 

To provide a comparison, we have performed calculations (using the NCP model) 
on the isostructural material H 3 [Co(CN) 6 ]. This is a useful comparison because 
experimental data for H3[Co(CN)g] show that that the thermal expansion is much 
smaller in this material (see table[l]). We used the same structure as for Ag 3 [Co(CN) 6 ], 
with the H lying exactly half way between the neighbouring N atoms. There has been 
some discussion in the literature concerning whether the positions of the H atoms 
might be disordered. In particular, the structure reported by Haser et al |25j has two 
sites for each H atom. However, our own recent neutron powder diffraction studies [2] 
have found no evidence for positional site disorder, with the best structure refinements 
having the H atoms in ordered positions, albeit with evidence for transverse vibrations 
as might be expected. 
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Table 4. Results for the calculated ground state of H3[Co(CN)6] and 
Au3[Co(CN)g]; atomic coordinates follow the definition in table [2 with H and Au 
replacing Ag, and experimental results for deuterated-H3 [Co(CNj6] were obtained 
at a temperature of 8 K [2]. 





d-Wo, [Co(CN)g] experiment 


H 3 [Co(CN) 6 ] calculation 


Au3[Co(CN)g] calculation 


a (A) 


6.412 


6.545 


7.990 


c(A) 


5.718 


5.678 


6.193 


C x 


0.238 


0.238 


0.196 


C z 


0.194 


0.198 


0.180 


N x 


0.376 


0.377 


0.312 


N z 


0.322 


0.324 


0.292 


C-N 


1.151 


1.161 


1.163 


H/Au-N 


1.291 


1.283 


1.978 


Co-C 


1.885 


1.917 


1.918 


H-H/Au-Au 


3.206 


3.272 


3.995 



The results of the DFT calculations are summarised in table [4] The agreement 
between the structures calculated by DFT and obtained from experiment is reasonable, 
and in line with the accuracy that is typical for DFT calculations: the errors on the 
values of a and c are +2% and —0.6%. In light of the above discussion on the role 
of Ag...Ag interactions, it is interesting to note this good agreement, which indirectly 
confirms the important role of Ag...Ag interactions in determining the ground state 
of Ag 3 [Co(CN) 6 ]. 

In figure [6] (left side) we show the contour plot for H3[Co(CN)g] analogous to 
the plot shown for Ag3[Co(CN)g] in figure [3] together with the locus of the minima 
points along the valley in the energy surface. In comparison with Ag3[Co(CN)g] we 
see a similar steep valley in the energy surface that provides the correlation between 
inverse changes in lattice parameters a and c. However, the changes in energy along 
the floor of the valley are significantly larger than in Ag3[Co(CN)g]. Given that this 
is the most significant difference between the calculations for the two materials, we 
surmise that it is closely associated with the differences in values of the coefficients of 
thermal expansion. 

3.7. Calculations on Au 3 [Co(CN) 6 ] 

Finally, we also report calculations of the ground-state structure of Au3[Co(CN)g], 
table [4] and a calculation of the energy surface is shown in figure [6j We believe that 
Au3[Co(CN)g] has yet to be synthesised, so these results are predictive. The energy 
surface is very similar to that of Ag3[Co(CN)g] (figure[3]), and thus we can predict that 
the thermal expansion will also show a similar correlation between a and c. Moreover, 
the variation of energy along the floor of the valley in the energy surface is very close in 
size to that of Ag3[Co(CN)g] - and is much smaller than that of H3[Co(CN)6] - so we 
anticipate that Au3[Co(CN)g] will large values of its coefficients of thermal expansion 
similar in size to those of Ag3[Co(CN)g]. Because there are no experimental data on 
this structure, it is impossible to judge the extent to which an additional dispersive 
Au...Au interaction is required to give a complete description of the energy surface, but 
the work of Pyykko and co-workers [20, ,2TJ|S2] suggest that there should be significant 
metallophilic interactions as in Ag3[Co(CN)g]. Furthermore, the calculated value of 
the a lattice parameter is much larger (of order 1.43 A) than the experimental value 
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Figure 6. Contour plot for the DFT calculations for Ji3[Co(CN)6] (top left) and 
Au3[Co(CN)6] (top right) with the trace of the energy of the locus of points along 
the bottom of the valley in the contour plot (bottom). It is interesting to compare 
these plots with the data presented in figure [3] 



for the isostructural La[Au(CN)2]3, where a — 6.664 A at 213 K [53]. This suggests 
that it is highly likely that the behaviour of Au3[Co(CN)g] and analogues such as 
La[Au(CN)2]3 will show colossal positive and negative thermal expansion. 



4. Analysis 

The DFT results shown in figure [3] point clearly to the existence of an energy surface 
in the space defined by the a and c lattice parameters that consists of a valley with 
very steep sides but only a shallow slope along the locus of minima points. The shape 
of the energy surface thus explains the origin of the negative thermal expansion along 
c for a positive expansion along a, consistent with the intuition from the structure 
outlined in the Introduction. 

We now consider the question of why the coefficients of thermal expansion of 
Ag3 [Co(CN)g] have colossal magnitudes. In our analysis we anticipate that the thermal 
expansion can be driven as a normal anharmonic phonon process with a very shallow 
energy surface. We write the free energy in standard form |27j as 

where V is the volume occupied by N atoms, Vq is the corresponding volume at the 
minimum-energy configuration, E is the potential energy of N atoms, and the phonon 
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free energy has been expressed in the high-temperature limit {k^T > TiuS) where u) 
represents an average phonon angular frequency. Minimisation with respect to volume 
yields 

V = V + 3NjBk B T = V (l + /3T) (2) 

where 7 is the normal overall Griineisen parameter, and the bulk modulus B = 
(Vod 2 E/dV 2 )~ 1 . The equation defines f3 = 3NBkBj /Vq. Our approach here is to 
use the energy calculations and experimental data for /3 (obtained as the trace of the 
linear thermal expansion coefficients) to obtain a value for Bj, and from an estimate 
for B (this can only be an estimate because of the use of the post hoc correction for 
the dispersive energy) to obtain an estimate of 7. Thus we obtain Bj = 47.3 x 1CP 12 
Pa -1 from the experimental data. The origin of the colossal values of the coefficients 
of thermal expansion can be expressed in terms of whether Ag3[Co(CN)g] has an 
unusually large value of B or 7 or both. 

From the energy calculations we obtained a fitted value for B^ 1 = 9.9 GPa, which 
would yield a value of 7 = 0.47. We stress that our fitted value of -B -1 is hard to 
estimate accurately from our DFT data because, being a small value, it is affected by 
very small errors in the calculations. However, it is nevertheless clear that the value 
for B~ x is small (by comparison, values of B^ 1 for NaCl, CaF2 and cordierite are 
24.4, 82.7 and 129 GPa respectively), so that B is large, and this had been anticipated 
from the significant strain broadening seen in diffraction experiments PQ. In fact, the 
value of B is so large that we have obtained a fairly typical estimate (in fact, on the 
low side of normal) for the value of 7 rather than a particularly large value. 

We note that we do not anticipate a large value for 7 for two reasons: first, 
because we do not anticipate the bond-stretching vibrational frequencies to change 
significantly with volume given that the bond lengths and bond orders do not change 
in the DFT calculation; and second, because negative contributions to the overall 
Griineisen parameter will also decrease its value. These will arise from dynamic bond- 
bending motions along the Co-C-N-Ag-N-C-Co chains, which will have the effect 
of pulling the end Co atoms closer together^]] We note that these motions have been 
observed in our RMC study. The final value of 7 will be a balance between motions 
that expand the volume and those that want to pull the structure inwards. 

What this analysis shows is that the colossal values of the coefficients of thermal 
expansion can be explained reasonably well with a normal, albeit slightly small, value 
of the phonon Griineisen parameter, and that the values of the coefficients of thermal 
expansion are so large because of the incredibly shallow energy surface along the 
bottom of the valley. 

We repeat this analysis for our results for H3[Co(CN)e]. We obtain B-f — 
8.36 x 10 -12 Pa -1 from experimental data. Our calculated energy surface gives 
B^ 1 = 56.4 GPa, from which we obtain 7 = 0.47. This value of 7 is very similar 
to that of Ag 3 [Co(CN) 6 ]. The principle difference between the thermal expansion of 
the two materials is thus due entirely to the different values of B, which is seen in the 
differences between the energy surfaces. 

% An analysis of the numbers of structural constraints and degrees of freedom shows that there 15 
more degrees of freedom than constraints per formula unit. This imbalance gives rise to rotational 
motions of stiff CoCg, CN and AgN2 units that will cause a shortening of the Co. ..Co chain. These 
motions correspond to the Rigid Unit Modes that are known to give negative contributions to the 
overall Gruneisen parameter [6] 1281 1291 130| . 
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Finally, from the shapes of the energy surfaces, and from the result that both 
Ag 3 [Co(CN) 6 ] and H 3 [Co(CN) 6 ] have normal values for 7, we make the prediction 
that Au3[Co(CN)g] will also show the same behaviour for its thermal expansion as 
Ag3[Co(CN)g], as also may other isostructural materials containing silver or gold. 
The main challenge involved in studying this compound experimentally is the limited 
aqueous chemistry of the Au + cation, which prevents a simple modification of the 
synthetic route for Ag 3 [Co(CN) 6 ]. 

5. Conclusion 

In this paper we have used DFT calculations to demonstrate that the coupled positive 
and negative thermal expansion in Ag3[Co(CN)@] and H3[Co(CN)g] arises from the 
existence of steep valley in the appropriate energy surface that can be traced to the 
existence of chains of atoms within the crystal structure that can flex like a garden 
trellis. The floor of the valley is particularly shallow in the case of Ag3[Co(CN)g], 
which gives rise to a large value of the volume compressibility (small value of the 
bulk modulus). This factor alone gives rise to the existence of colossal values of 
the coefficients of thermal expansion, with the values of the Griineisen parameters 
appearing to be normal. This work has pointed to the important role of Ag...Ag 
metallophilic interactions, although we cannot claim to have arrived at a complete 
understanding of how these give rise to such a shallow floor of the energy surface. 
More work is needed in this regard. Moreover, we have predicted that Au3[Co(CN)g] 
will show similar behaviour to Ag3[Co(CN)e], and we anticipate that other materials 
with metallophilic interactions will also display extreme behaviour such as collosal 
thermal expansion. 
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